FINITE ELEMENT ANALYSIS OF MAGNETIC DAMPING EFFECTS ON G- 

JITTER INDUCED FLUID FLOW 

1 

Bo Pan and B. Q. Li ^ , 

Department of Mechanical Eng inee ring 
Louisiana State University 
Baton Rouge, LA 70803 

and 

H. C. de Groh 

NASA Lewis Research Center 
Cleveland, OH 44135 


ABSTRACT 


This paper reports some interim results on numerical modeling and analyses of magnetic 
damping of g-jitter driven fluid flow in microgravity. A finite element model is developed to 
represent the fluid flow, thermal and solute transport phenomena in a 2-D cavity under g-jitter 
conditions with and without an applied magnetic field. The numerical model is checked by 
comparing with analytical solutions obtained for a simple parallel plate channel flow driven by g- 
jitter m a transverse magnetic field. The model is then applied to study the effect of steady state 
g-jitter induced oscillation and on the solute redistribution in the liquid that bears direct 
relevance to the Bridgman-Stockbarger single crystal growth processes. A selection of computed 
results is presented and the results indicate that an applied magnetic field can effectively damp 
the velocity caused by g-jitter and help to reduce the time variation of solute redistribution. 


1. INTRODUCTION 

It is well known that natural convection affects the diffusional growth of crystals from 
melts. Microgravity environment in space vehicles has provided a unique means by which the 
unwanted convective flows induced by gravity can be significantly reduced. However, the space 
environment often reveals effects that are not known or present on Earth. One of these effects is 
g-jitter or residual acceleration phenomena associated with microgravity environment. Residual 
accelerations or g-jitter occurring during space processing of melts, such as the melts from which 
sn^!e crystals are grown, can cause appreciable convective flows, making it difficult to realize 
diffusion controlled single crystal growth in space[l]. 

G-jitter comes from crew motions, mechanical vibrations (pumps, motors, excitations of 
natural frequencies of spacecraft structures), spacecraft maneuvers and attitude, atmosphere drag 
an the Earths gravity gradient [2], Studies on g-jitter effects indicate that convection in 
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microgravity is related to the magnitude and frequency of g-jitter and to the alignment of the 
gravity field with respect to the growth direction or the direction of the temperature gradient [3- 
6]. A residual gravity of 10' 5 where & is the Earth gravity constant, is sufficient to cause 
unacceptable fluid motion in the liquid, unacceptable in the sense that the g-jitter induced flow is 
intensive enough to affect deleteriously the solute distribution in the melt and hence the quality 
of the crystal grown. The orientation of the gravity vector with respect to the temperature 
gradient plays also an important role in melt flows. The velocity attains a maximum when the 1 
gravity vector is perpendicular to the temperature gradient. J 

Numerous attempts have been made to estimate and calculate the adverse effects of time 
varying g-jitter [6-1 1]. These calculations are all based on numerical models. Both 2D and 3D 
numerical models have been developed for this purpose [4,12]. These models have been used to 
study the effects associated with both idealized single- and multiple-frequency g-jitter 
modulations and realistic g-jitter data collected by accelerometer during actual flight 
experiments. These studies showed that the frequency, amplitude and spatial orientation of the 
residual gravity vector all play an important role in determining the convective flow behavior of 
the system. When the residual accelerations oscillate about the positive and negative of an axis, 
the orientation of this direction relative to the density gradient determines whether a mean flow is 
generated in the system [11]. Sinusoidal oscillating accelerations induce an oscillating 
convective flow and composition oscillation in the liquid. The resultant velocity and 
composition fields oscillate with the same frequency of the affecting gravity field. 


These studies have shed light on the basic nature of g-jitter effects and have provided a 
thrust to devise useful mechanisms to suppress g-jitter induced convective flows. One technique 
is to apply the principle of magnetic damping to achieve the velocity reduction. Magnetic 
damping is less effective than reducing gravity for convection reduction. It originates from the 
interaction of the flow field of an electrically conducting fluid and an applied magnetic field. 
This interaction gives rise to an opposing Lorentz force that reduces or damps the flow velocities 
in the liquid. Application of a magnetic field has been found to be effective in controlling 
convection during crystal growth from melts under terrestrial conditions and has now been 
widely practiced in the metals and semiconductor industries [13]. While there has been a large 
amount of literature addressing both the basic characteristics of g-jitter induced flows and 
application of magnetic damping to convection control in melt growth on Earth, little research 
work has been done on the assessment of .the magnetic field effect on the convective flow 
associated with g-jitter in microgravity. 


In this paper, a numerical analysis is presented of g-jitter driven flows and their effects on 
the solutal element distributions in a simplified Bridgman-Stockbarger crystal growth system. 
The numerical model is based on the solution of the Navier-Stokes equations for fluid flow, the 
energy balance equation for the temperature distribution and the mass balance equation for the 
species transport in a Ga-doped germanium melt. The model is applied to study the magnetic 
damping effects on the oscillating convective flows driven by g-jitter and on the solutal 
redistribution in the melts. Both single and multi-frequency g-jitter effects are studied and 

numerically obtained solutions are compared with available analytical solutions whenever 
possible. 
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2 . PROBLEM STATEMENT 

^Vr temperature T h . The lower side walls are adiabatic awT^^JSj! 1 . * 

constant velocity V m at the top of the computational space (inlet) The amnnnl * f - Pm a a 
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Figure 1. Schematic representation of the natural convection 
problems in a 2-D cavity in microgravity environment 

is environment, the residual acceleration experienced by the spacecraft 

time ST a °5 3 Cady com P° nent ^ a time-dependent component. The magnitude of the 

conside^e ^eT^TV 8 T* ?""** ^ ** ° ffhe ^ component. In Ais pape^ we 
time dependent acceleration (g-jifter) represented by Ae following expression: 

SW = Z sin(2ar /„/) e .r '> * 1 e< l n 

n-i 

res Wuaf gra^tv 1 is^ara &e ? Uency of g-J itter respectively. The direction of Ae 

77 terms w * 5**“* interf f ce > ^ has Ae most detrimental effects in J 

^ ~magiiitkl^<r^ ch is T nf ^J ab ° V ? CITStal ^^dilystem is immersed m“a transverse 
magnetic held, which is intended to damp Ae convection resulting from g-jitter. 

field onTS in b th CtlVe ^ !° Stady u g ' jitter ' driven melt flow ^ Ae effect of Ae magnetic 
““ . u AC CfyStal gr0Wth system ^ described above. To simplify Ae analysis Ae 
sojidification phenomena are not considered and Ae melt-crystal in terface is heldfl gr The melt 

^ <7^j- ^v,y, 
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fnn'rinS 2 '° 35 ilIustrated in Figure 1 are governed by the 

ty equa ion, the Navier-Stokes equations, the energy balance equation and the species 

folto™ Cq ° n ' 111686 equations ma y be ^tten in a nondimensionalized form as 
Equation of_ Continuity- 


V»u = 0 

Equation of_ Momentum Balance: 


( 2 ) 


du 


— + u • Vu = -Vp + PrV 2 u - RaPrg*T + Ha 2 Pr(u x b x b) 


(3) 


Equation of Energy Balance 


dT ? 

— +u*vr = v 2 r 

at 


(4) 


Equation of Solute Transport- 


Scf 3C 

Prl^/ 


+ u • VC 


= v 2 c 


(5) 


where u-O^ig is velocity field, p the pressure, b applied magnetic field, T the temperature, C 
the concentration of solutal element. The last term in Equation (3) represents the effect of the 
Lorentz force that arises from the application of an external magnetic field. In the above 
equations, the velocity is scaled by k/ L, the time by L 2 /k, the pressure by P k 2 /L 2 . The 
temperature is defined by (T(r,t)-T m )l(T h -T m ) where T(r,t) is the dimensional temperature. 
The dimensionless parameters, that is, Ra the Rayleigh number, Ha the Hartmann number, Pr 
the Prandtl number, Pe the Peclet number and 5c the Schmidt number, are defined as: 


Ra 


Ha = B 0 L 





V 


D ’ 


The boundary conditions for the system under consideration are given by: 


T = 1, 




Pe g pr 

Sc ’ 


uxn = 0, 


dC 

— Peg(l — k)C; at z=0 


T = 0, 


u • n 


Pe g Pr 


u x n = 0, 


dC 

~^=Pe g (C-l); at z=l 


5c 
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u*n = 


F ‘x Pr 


Sc * u * e M' e w .VC-0; along the two sides 
r = l along the hot zone side; and e„.vr =0 along the adiabatic zone 

where n is the unit vector normal to the planar crvstal-melt inww 

wall, and k is the solute segregation coefficient ' ’ " 0nnal t0 1116 side 

3. FINITE ELEMENT SOLUTION 

Galerkin finite element melh^T^stifeess^tok i^b^dT^^ 003 ^ SOlVed US “ 18 ±C 
Weighted Residuals. Within each element die by mmg Gderkin ’ s method of 

shape functions of* vi/ a^d 0 TTienA?^ VanabIes u ’ P > T > C ** interpolated by 

are approximated by- * " ^ prCSSUre ’ and concentrTon fields 


u i(x,t) = <t> T \j .(/) 
F(x,r) = ^P(,) 
F(*.0 = 0 r T(O 
C{x,t) = ^ C (/) 


( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 


where U„ P, T and C are column vectors of element nodal point unknowns. 

R anH S R tUtm |?- ^ ab ° Ve eqUations into governing equations, we get the residuals of R R 

errors to zero and the shane fimrtir.no u , ghted Residuals seeks to reduce these 
Ute procedure's given in [1 7), the governing'.^ F °" OWi ” g 

Momentum 

( £ <M?dV) ~j- + ( gtu .Vt T dV)V, - ( J[ itr Vj T dV)P 

+ ( i Pr V?s • r ^)U/ + ( ( PrV,< . U/ + ( J[ fe Pr 

= £ | /fa 2 Pr(u x b x b)rfF 

Conti nuity ( penalty fnrmi.latmn ) 

( J[ y/V <f> T dV)Vj =-e(j v v/y/ T dV)P 

Energy 


( 10 ) 


( 11 ) 


([M T dV) ^ +(l0u»V0 r dV)T + (j v V0.V0 T dV)T = - facets 


( 12 ) 
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Species 

( l P^ eeTdV> > + C *' + ( j[ V0* V6 T dV)C = - [qjds 


( 13 ) 


Once the form of shape functions <J>, 0, and y are specified, the integrals defined in above 
equations can be expressed by the matrix equation. Combining the momentum and energy 
equations mto a single matrix equation, we get the element stiffoess matrix equation: 


M 0 0 Tuj rA(U) + K + jEM;'E r B 0 Tu 

0 N r 0 It + 0 D r (U)+L r 0 T 

.0 0 N cJ_cJ [ 0 0 D c (U) + L c ]_C 


F 

G r 

9c 


(14) 


Note that in constructing the above element matrix equation, the penalty formulation has been 
applied, and P in the momentum equation is substituted by|M;'E r U. The assembled global 

matrix equations are stored in the skyline form and solved using the Gaussian elimination 
method. The coefficient matrices of above are defined by: 


j* 

II 

s 


M p = \ v yny T dV 


N r = | ee T dv 

N c = [y T ee T dv 

K= j[Pr V<j>vfdV 


E = \ v {V<t>)y/ T dV 


L r = j V6*V0 T dV 

L c = \ v9*VG T dV 

A(U)= [tvVfdV 


D r (U)= \ v evve T dv 

V c (\})=\^eu.Vd T dV 

B = ^ RdPxg<f>0 T dV 


F = £ r<f>dS + Ha 2 Pr(u x b x b )dV 


Gj = -^qjOdS 

Gc = — J q c 6dS 


4. RESULTS AND DISCUSSION 

The finite element model described above enables the prediction of the transient convective 
flows driven by g-jitter, the temperature distribution and the solutal redistribution in the melt 
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with and/or without an applied magnetic field. The calculations were performed using a finite 

(oSiritl'th^T deVe ‘T d by 0Ur reSeaiCh The code employs the backward Eder 

(or implicit) timescheme, with automatic time step control, and the penalty method for pressure 

r mP , UtCr C ,° de ™ Ch6Cked against other mailable commercial flow codes, 
, . 8 DAP > analytical solutions for a wide variety of flows and thermal problems A 

327 rr, ed w " r ^ * e parameters — * t 

Tab \ T? e traxisicnt calculations are started with an initial thermal and fluid flow field 
obtamed with a constant gravity field; for the results in Figures 3-10 the initial flow field is 

uTth^crystal-mdt interface.^ ^ “ H, “" 11 “ d ‘2 -e sorted wid, g =10 ^, both parallel 



Figure 2. Comparison of the numerical and the analytical solutions for 
convective flows in a parallel plate channel: Y the location between the 
plate and U the fluid flow velocity, both nondimensionalized. 


fl T °, ^7^ ass ® ss 4116 code P redlcta bility for the magnetic damping effects on g-jitter driven 
flows the finite element results are compared with the analytical solutions obtained for 
convec ive flows in a simple parallel plate geometry subject to a transverse magnetic field [1 81. 
The height-to- width aspect ratio is set at 5, which was found to be sufficient to simulate the 
infinitely long geometry. The results are shown in Figure 2. Clearly the numerical results 
compare veiy well with the analytical solutions. 




Table 1; Pa rameters used in calculations m 
Property 


Thermal conductivity of the melt (k) 
Heat capacity of the melt (Cp) 
Density of the melt (pj 
Kinematic viscosity of the melt (v) 
Solute diffusivity (D) 

Thermal diffusivity of the melt (k) 
Thermal expansion coefficient (P) 


17 WW 

0.39 J g' 1 K' 1 

5.6 xlO 3 Kg m’ 3 

1.3xl0' 7 m 2 s' 1 

1.3xl0 4 mV 

1.3xl0* 5 m J s* 1 

2.5xlO‘ 4 K*' 


Value 


Prandtl number (Pr) 

<X0l 

Peclet number (Pe) 

5.0 

Rayleigh number (Ra) 

1.45xl0 3 ( on earth value ) 

Hartmann number (Ha) 

10-100 

Schmidt number (Sc) 

10 


Operating conditions Value " 

Hot zone temperature (T h ) 133 1 K 

Distance between inlet and interface (L) 0 . 01 m 

Interface temperature (T.J 1 23 1 K 

Translation rate (VJ 6.5 xlO 4 m s' 1 

Gravity (g)/ Frequency (Hz): 1 O' 2 , 1 0\ 10^1 ,0.1, 0.01 

Magnetic field B 0 0.022-0.22 T 

* lgo = 9.8 ms' 2 



Figure 3. Temperature distribution in the 2*D cavity; thermal transport in the 
system is primarily attributed to conduction. 
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Figures 4 and 5 show the basic fluid flow structure development in the melt pool with and 
without an applied magnetic field. The magnetic field effect on the g-jitter induced flows is 
er i ustrated in Figure 6, where the maximum velocity (U component) in the pool is plotted 
as a function of time for the cases with and without an applied magnetic field. For the case 
consi ered, an application of a somewhat moderate magnetic field results in a reduction of 
ve ocity level by almost 50%. The magnetic damping effects are also related to the applied 
magnetic field strengths and the field directions, as shown in Figure 7. The higher the applied 
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during crystal Growth* 1 To f'^ U ^.^°°* ^* ve a significant effect on the solutal distribution 

SSttZ “ mves " g f : *“ effea - «* *">« dependent concentration profiles in 
liquid pool are plotted tn Figures 8 and 9. From Figure 8, it is seen that the concentration 
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evolvesastime continues and the g-jitter induced flow strongly affects the solutal transport in the 
melt. With an applied magnetic field, the liquid convection is reduced and hence its effect on the 
solutal distribution is also reduced, as appears in Figure 10. Figure 10 compares the time 
variation of the concentrations at two points in the melt pool. Apparently, with an applied 
magnetic field, the concentration oscillation is also greatly reduced . 



(a) (b) 

Figure 6. Variation of maximum flow velocity in the Ga-doped germanium melt as affected by g- 

jitter and an applied magnetic field: (a) g=10 J , 10°, lfr 4 , fii= 1,0. 1,0.01, B=0; and (b) g=10 1 , 10 \ 

10 -4 , fii= 1,0. 1,0.01, Bx=0.22 Tesla, Ha=100. Gravity oscillates parallel to the crystal-melt 
interface. 



(a) (b) 

Figure 7. Dependency of transient maximum flow velocity in the Ga-doped germanium melt on the 
direction and magnitude of the applied magnetic field: (a) effect of applied field strengths and (b) the 
effect of magnetic field directions. Other conditions: single frequency K-iitter e=10' 3 and f =0 1 
acting in the x-direction. 
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5' e . c u rystal . gr0Wth P 0 '"' of ™ w . *e concentration distribution along the 
R^l ^l anuoportant implication to the impurity striation in the crystals gLvn. 

Sti aTone L' ^T V a 6 ^ ° f ** app,icd “*"* field on the concenLion 
field fst Figme m Z ^ Wi,h the ab “"“ of an external magnetic 

tele ^ 4^ HifT ^ “" tratl °“ pr0file osciUates wi* time along the interface and as 

i ndla,;„p thet^nTeff* ?T‘‘ Cmin/Cm ' m,) “ ““““rati 0 " is observed at some instances, 
indicating the strong effect of convection m the liquid pool. With an applied magnetic field, 
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however, the oscillation is peatly reduced, and the concentration profile along the interface 
evolves over a much longer time period, suggesting that the convection effects are much reduced 
Comparison of Figures 11 and 12 shows also that the oscillation in concentration distribution 
along the interface is reduced as a result of the application of a magnetic field. 



Figure 9. Concentration profile of multi-frequency Figure 10. History of concentration 

g-jitter, with magnetic field Bx=0.22, Ha=100. change with time with and without 

The profile remains approximately the same for magnetic field at certain points, 

the time period as shown in Figure 6. 



X Axis 


Figure 1 1 . Concentration distribution along the 
interface at different times without imposed 
magnetic field. 



X Axis 


Figure 12. Concentration distribution along 
the interface at different times with a 
magnetic field Bx=0.22, Ha=100. 






5. CONCULDING REMARKS 

growth processes. The model^predictions compared^we^S^I 1 ^ 1 ''’ 3 ' 8 ^ 5 ' 1181 ' Cr ?' stal 
The numerical model was applied to studv the ft m n . . avai a ^ e analytical solutions, 
solute redistribution. Z ZtericL r^ts obta^ r ^ h “‘ <"«* and 

magnet, c field can have aTJfilt effec^oTth model sho ^ that an applied 

redistribution in the system. C convectIve R° w and hence the solute 
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